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One of the most experimentally testable explanations for the origin of the baryon asymmetry of the 
universe is that it was created during the electroweak phase transition, in the minimal supersym- 
metric standard model. Previous efforts have focused on the current for the difference of the two 
Higgsino fields, Hi — H2, as the source of biasing sphalerons to create the baryon asymmetry. We 
point out that the current for the orthogonal linear combination, H1 + H2, is larger by several orders 
of magnitude. Although this increases the efficiency of electroweak baryogenesis, we nevertheless 
find that large CP- violating angles > 0.15 are required to get a large enough baryon asymmetry. 
PACS: 98.80.Cq 



1. A highly constrained proposal for the origin of bary- 
onic matter in the universe is electroweak baryogenesis 
in the minimal supersymmetric standard model (MSSM) 
• Unlike many other baryogenesis mechanisms, this 
one has strong prospects for being falsified in upcoming 
experiments, due to its need for some light exotic parti- 
cles, notably a right-handed top squark which is lighter 
than the top quark fj]J|. 

The basic mechanism Q is intuitively clear: parti- 
cles interact in a CP-violating manner with bubble walls, 
which form during the first order electroweak phase tran- 
sition, when the temperature of the universe was near 
T = 100 GcV. This causes a buildup of a left-handed 
quark density in excess of that of the corresponding anti- 
quarks, and an equal and opposite right-handed asym- 
metry, so that there is initially no net baryon num- 
ber. The left-handed quark asymmetry biases anoma- 
lous sphaleron interactions, present within the standard 
model, to violate baryon number preferentially to create 
a net quark density. The resulting baryon asymmetry of 
the universe (BAU) soon falls inside the interiors of the 
expanding bubbles, where the sphaleron interactions are 
shut off, and thus baryon number is safe from subsequent 
sphaleron-induced relaxation to zero. 

Despite the simplicity of this idea, a quantitatively ac- 
curate calculation is difficult to carry out. Various ap- 
proximations have been made, leading to a variety of 
formalisms which give somewhat conflicting predictions 
for the dependence of the BAU on the parameters of the 
MSSM. Although most authors agree on the diffusion 
equations which govern the generation of the left-handed 
quark asymmetry, there is less consensus about how to 
derive the source terms which appear in these equations. 

In a previous paper 0, we advocated an approach 
based on the classical force on particles fl(]Jll| due to 
their spatially varying masses as inside the bubble wall. 
It is straightforward to solve for this force, put it into 
the Boltzmann equations, and derive corresponding dif- 
fusion equations. No ad hoc assumptions are needed; one 
only requires that the width of the wall be significantly 



larger than thermal de Broglie wavelengths of particles 
in the plasma, to justify an expansion in derivatives of 
the background Higgs field, which constitutes the bubble 
wall. Detailed calculations of the wall width confirm that 
this is indeed a good approximation p"2[ . 

One puzzling conflict between our earlier work H and 
that of others was that we derived a source which remains 
large when the ratio of the two Higgs fields, H2(x) /H\ (x 
is constant inside the bubble wall. Other authors J 
found sources proportional to the derivative of this quan- 
tity. Careful analyses of the shape of the wall have shown 
that in fact H2/H1 remains constant to within a part in 
10 2 - 10 3 fHH , so that the dependence on d(H 2 / 'H x ) / 'dx, 
if correct, would result in a large suppression of the gen- 
erated BAU. 

In this Letter we explain the origin of the apparent 
discrepancy: it is the result of different choices about 
which linear combination of Higgsino currents is used to 
source the diffusion equations. References consid- 
ered only the combination Hi — Hi , based on the obser- 
vation H that the orthogonal source, H± + H2, is driven 
to zero in the limit that interaction rate T proportional 
to the top quark Yukawa couplings becomes infinite. Al- 
though this approximation is convenient because it sim- 
plifies the network of coupled diffusion equations, in the 
present application it can lead to a serious underesti- 
mate of the BAU. Parametrically, the Hi + H2 source 



is suppressed by {FDjJ 



-1/2 



20/T is the 
However for realistic 



where 



Higgsino diffusion constant 14 1. 

values, (TDh)^ 1 / 2 ~ 1, so there is no actual suppres- 
sion, eliminating the need for large CP violating phases 
to get the observed BAU, which we had found with the 
Hi — H2 source j|. Here we will demonstrate how this 
comes about, while at the same time updating our ear- 
lier computation so as to give a quantitatively accurate 
determination of the BAU as a function of relevant pa- 
rameters in the MSSM. 

2. To explain the how the classical force mechanism 
works, we first review the simple case of a top quark 
with a spatially varying complex mass [fill, m(x) = 



1 



yH(x)e l8 ( x \ By solving the Dirac equation (i$ — 
\m\ cos9 — i\m\ sin #75)?/; to first order in derivatives (the 
WKB approximation) , one finds that a particle of energy 
E experiences a spin (s) dependent force 
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The spin dependent part of the force has opposite sign 
for antiparticles, which causes the distributions of likc- 
helicity fermions and antifermions to be separated in the 
vicinity of the wall. For relativistic particles, we can 
approximate helicity by chirality, and speak of spatially 
varying chemical potentials, fJ. L ,R for the left- and right- 
handed components; these are related to the asymme- 
try between the particle and antiparticle densities by 
n(x) — n(x) = /i(x)T 2 /6. Diffusion equations can be de- 
rived by inserting the force (jj]) into the Boltzmann equa- 
tion, doing an expansion in moments of the distribution 
functions, and truncating the expansion. They have the 
form 



Dfi L 



v w (J, L 'i + ^^Li = Si(x), 



(2) 



where D is the diffusion coefficient (of order the inverse 
mean free path), v w is the bubble wall velocity, and T is 
a damping rate representing decays or inelastic collisions 
of the left-handed fermions. The source term for each 
species is related to the classical force through a thermal 
average H 



S(x) 



v w D 



(v x F(x))'. 



(3) 



Although this looks similar to spontaneous baryogenesis 
jL5|, it is not the same; the latter arises through the col- 
lision term in the Boltzmann equation, while ours comes 
from the flow term [||JlJ], and we find that it gives a 
somewhat larger effect. 

Once the left-handed quark density is found from eq. 
(0) , the density of baryons generated by sphalerons can 
be computed by integrating [i L (%) = E gi M^i in front of 
the wall: 
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where r sp h is the diffusion rate of the Chcrn- Simons 
number as measured by lattice simulations Jl6| : r sp h = 
(20 ± 2)cq v T. The exponential accounts for sphaleron- 
induced relaxation of the baryon asymmetry back to zero 
due to restoration of thermal equilibrium, in the limit of 
a very slowly moving wall. The coefficient depends on 
the squark spectrum: it equals to 45/4 if all squarks are 
heavy and 72/7 if only the right handed stop is light. 

3. When the above picture is adapted to the MSSM 
some complications occur, because the quarks do not get 
complex masses and hence are not directly sourced. The 



chargino mass matrix does contain complex phases how- 
ever, and the chiral asymmetry which develops in the 
chargino sector induces one for the quarks because of the 
strong coupling between Higgsinos and the top quark. 
This system involves not just one, but many coupled 
diffusion equations. We will show that they can, nev- 
ertheless, be reduced to a single equation by reasonable 
approximations. 

We start the treatment of the MSSM by deriving the 
source term analogous to S(x) in eqs. (§-^)- The mass 
term in the Lagrangian for the charginos is 



i> R M x i) L = (uS 
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plus the Hermitian conjugate, ip L M^ipR. Because there 
are two Higgs doublets, Hi and H 2 , there are two corre- 
sponding Higgsino fields h 12 - w is the wino, superpart- 
ner of the W boson, and g is the weak gauge coupling. 
The complex phases of the wino mass m 2 and the fi pa- 
rameter are the origin of a CP violating force. By again 
solving the Dirac equation in the WKB approximation, 
one can find the spin-dependent forces which act on the 
two mass eigenstates, with masses m\. Since Higgsinos 
couple strongly to top quarks, we are interested in the 
one which smoothly evolves into the Higgsino state in 
front of the wall, where -Hi, 2 —* 0. The spin-dependent 
part of the force has the same form as in eq. ([!]), with 
the spatially varying phase given by 



g 2 lm(m 2 [i) 
(m 2 + — m 2 _) 



{H\H' 2 + H[H 2 ) 



(6) 



where the Higgsino-like mass eigenvalue, m~(x), is m 2 _ 
(the lighter one) if |//| 2 < |m 2 | 2 , and m 2 + otherwise. 
We remark that the effective WKB expansion parameter 
O'h/E remains small even when mass gap m 2 ^ — m 2 _ at- 
tains its minimum value; this corresponds to \m 2 \ — \fi\, 
where one has parametrically 0' h /E ~ gHi/wE\fi\ <C 1, 
because for typical wall widths wE > 20 fl2] ]. 

Having specified the source term, we must now deal 
with the diffusion equations in which it appears. In gen- 
eral, these are a set of coupled equations for the two Hig- 
gsinos (ft.1,2), the winos, 6 flavors of right-handed quarks, 
3 generations of left-handed quark doublets, and all the 
superpartners which are light enough to be present at 
the temperature T. Following Huet and Nelson S, we 
will make the simplifying assumption that the super- 
gauge interactions {e.g., the coupling of winos to quarks 
and squarks) are in thermal equilibrium, so that particle 
and corresponding sparticle chemical potentials are equal 
to each other, species by species; in particular gaugino 
chemical potentials are driven to zero. We will further 
assume that the right-handed top squark is the only light 
squark. The latter must be light to get a strong enough 
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electroweak phase transition, and it turns out that mak- 
ing all the others heavy is favorable to electroweak baryo- 
genesis, as well as satisfying constraints on CP violation 
in the MSSM. 

With these simplifications, the diffusion equation net- 
work can be reduced to four equations, for the poten- 
tials of the Higgsinos hi, h 2 and the left-handed third 
generation quark doublet and the right-handed quark 
t#. Ignoring interactions involving left-handed squarks, 
which we assume to be heavy, the important interactions 
coupling these species come from the potential 



V = yq 3 H 2 t R + fih\h 2 + h.c. 



(7) 



The rate for the Yukawa interaction is T y . In addition 
there are helicity flipping interactions with rate Th f cou- 
pling h\ and h 2 , due to the /i term, and in the broken 
phase inside the bubble (x < 0), the top quark Yukawa 
coupling becomes the top mass, which causes helicity 
flips between the left- and right-handed top quarks, q% 
and tn, with a rate of T m . Defining the diffusion op- 
erator T>i = —6(Did% + v w d x ), and defining also the 
convenient linear combinations /i± = ^(/J^ i ) an d 
H V =H + - /i_ - /if „ + /i 93 , we have 

V h fi + + 2T hf fi + + ^Tylly = S H (8) 

2? h M- - \T v iiy = (9) 

\T> q \x qz + ^n~ h Ty^y - f m + 2T SS = (10) 

\V q n tR - ^n^Tyfiy + f m - r ss = 0. (11) 

In deriving these equations we have assumed that all 
squarks but right handed decouple from thermal equi- 
librium. The factor ~ ^x-Kiixj^, where x~ h = m^/T 
and K 2 is the modified Bessel function, accounts for par- 
tial decoupling of higgsinos. The top-quark helicity-flip 
term is defined to be T m = T m 0(—x)(^t R — Hq 3 )/T and 
that for the strong sphaleron rate is T ss = r ss (20// q3 + 
26fi tR )/T with T ss = 1500K sph a^T @. 

It is noteworthy that only the linear combination /z+ 
gets directly sourced in our WKB treatment, whereas the 
source for the combination /i_ vanishes (in this respect 
ref. JB| was in error). This is in contrast to papers which 
treat the particle reflections from the wall quantum me- 
chanically; these works find that both /x + and are 
sourced. Here we point out that, within the quantum 
reflection formalisms, it is still true that fi+ has a larger 
source than //_, because H[H2 + H' 2 H\ is larger than 
H[H 2 — H' 2 Hi fll8|| . However nobody has heretofore con- 
sidered the former because of the unfortunate approxima- 
tion of imposing equilibrium of the Yukawa interactions 
(T y — > oo), which forces fi + to zero. 

The network of diffusion equations can be integrated 
numerically to to obtain /i tR and fi q3 . From these one can 
compute the left-handed quark potential [II = 18/i* R + 
15/i 93 , which must then be numerically integrated in eq. 
(0) to obtain the baryon asymmetry. 



4. We have carried out the above procedure to study 
how the BAU depends on the velocity of the bubble wall 
v w , the wall width w (appearing in the Higgs field pro- 
file as i? (l - tanh(a:/uO)/2), tan/3 = (H 2 )/(Hi) (at 
zero temperature) and the chargino mass parameters fj, 
and m 2 . We take as our preferred values tan/3 = 3, 
w = 6/T Jl|, v w = 0.1 and fi = m 2 = 100 GeV. 
We also took T = 90 GeV and H = 110 GeV. We 
computed the values of D^, T y , T^f and T m from 
the actual Feynman diagrams and thermally averaged 
cross sections, and also using the results of ref. [ ^9| . 
With our preferred parameters we find Dr — 80/T, 
r m = 0.0012T, T h f = 0.001a:?/[(l + x h )(2 + x~ h )] and 
T y = [0.012 + 0.022 a;,, (1 - 1.87/4 + 0.05/4) 3 / 2 ]T. In 
the following figures where certain of these quantities are 
varied, those which are not specified have the above val- 
ues. It is customary to express the BAU as a scaled ratio 
of baryons to photons, 7710 = 10 w nB/n 1 = 7 x 10 10 ns/s, 
where s = (2n 2 /45)g*T 3 , and we take the number of de- 
grees of freedom in the MSSM to be g* = 110. Current 
limits from big bang nucleosynthesis give 3 7710 ^ 4. 

In figure 1, the BAU is plotted as a function of wall 
velocity, varying the width of the bubble wall to obtain 
the different curves. We have taken the CP-violating 
phase appearing in lm(m 2 jj,) (eq. @) to be maximal; to 
satisfy r/io = 3, one should rescale this phase accordingly. 
The efficiency of baryogenesis tends to peak for v w — 
0.01 — 0.1. Interestingly, recent work on bubble expansion 
has suggested just such a range of small values in the 
MSSM |20|. 




Figure 1: Baryon asymmetry 7710 as a function of wall velocity 
v w for wall width w = 6/T, 10/T and 15/T. 

We can also consider the dependence of the BAU on 
the chargino mass parameters, /i and m 2 . Figure 2 ex- 
hibits contours of constant 7710 assuming the maximal 
magnitude of the CP-violating phase, &rg(m 2 (i) = tt/2, 
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and wall velocity v w = 0.1. Since we need 7710 = 3, and 
the maximum value in the figure occurs around 7710 = 20 
for /i ~ iri2 ~ 50 GeV, we find that CP phase can be 
no smaller than 0.15, and then only if the charginos are 
rather light. Such large values of the CP phase could 
come into conflict with the most stringent present ex- 
perimental constraints, which come from searches for the 
electric dipole moment of the 199 Hg atom pl[ . In order 
to evade these constraints, one must assume that most 
of the squarks or gauginos contributing to EDM loop di- 
agrams are heavy, in the TeV range. 
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Figure 2: Contours of constant 7710, assuming maximal CP 
violation (arg(m2^) = 7r/2) and v w — 0.1, in the plane of the 
chargino mass parameters \x and m,2 (in GeV units). 

5. In summary, we have presented a quantitatively ac- 
curate analysis of the baryon asymmetry of the universe 
in the MSSM, using the classical force mechanism, which 
allows a consistent computation of the sources appear- 
ing in the diffusion equations. A CP-violating force acts 
on Higgsinos while they cross the bubble walls causing a 
particle-antiparticle separation in charginos, which then 
gets partially transformed to a chiral quark asymmetry 
that biases sphalerons to produce baryons. Unlike previ- 
ous authors, we have not assumed that the two compo- 
nents of the Higgsino, h\ and h^, reach complete chemical 
equilibrium through Yukawa and helicity-flipping inter- 
actions, and we showed that a large enhancement of the 
BAU can result, relative to what one would get this over- 
simplification. Nevertheless, after carefully computing 
the relevant damping rates and diffusion coefficients, and 
accurately solving the transport equations, we find that 
it is difficult to generate the observed baryon asymmetry 
unless CP violation in the /j, term is close to maximal. 
Further details of our formalism and calculations can be 



found in ref. |Q. 
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